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Abstract 

We consider the minimal supersymmetric triplet seesaw model as the origin of neutrino masses 
and mixing as well as of the baryon asymmetry of the Universe, which is generated through 
soft leptogenesis employing a CP violating phase and a resonant behavior in the supersymmetry 
breaking sector. We calculate the full gauge-annihilation cross section for the Higgs triplets, 
including all relevant supersymmetric intermediate and final states, as well as coannihilations 
with the fermionic superpartners of the triplets. We find that these gauge annihilation processes 
strongly suppress the resulting lepton asymmetry. As a consequence of this, successful leptogenesis 
can occur only for a triplet mass at the TeV scale, where the contribution of soft supersymmetry 
breaking terms enhances the CP and lepton asymmetry. This opens up an interesting opportunity 
for testing the model in future colliders. 

PACS numbers: 98.80. Cq,12.60.Cn,12.60.Jv 



Leptogenesis is an elegant way to generate the baryon asymmetry of the Universe in 
connection with the origin of the observed neutrino masses and mixing through the seesaw 
mechanism Q|. One way of understanding a tiny neutrino mass is to relate it with the small 
vacuum expectation value of a Higgs triplet |2j whose decay can also induce the cosmological 



mryon asymmetry in the presence of at least two Higgs triplets [3] or a right-handed neutrino 
^| as required by the generation of non-trivial CP and lepton asymmetry. In the minimal 
super symmetric version with one pair of triplets, there is a new way of leptogenesis (called 
"soft leptogenesis") in which CP phases in the soft terms can contribute to generate the 
.epton asymmetry HQ. Soft leptogenes, . in the minima. S up_et r ie Higgs tnp.et 
model has been considered first in Ref. [7(. 

In this paper, we revisit this last scenario to provide a careful analysis on the quantities 
for the lepton and CP asymmetries and their cosmological evolution by considering the full 
set of Boltzmann equations including thermal masses and the temperature supersymmetry 
breaking effects consistently. We will also derive a set of simple Boltzmann equations from 
the Maxwell-Boltzmann approximation taking into account the difference between the Bose- 
Einstein and Fermi-Dirac statistics, and show that they provide a fairly good approximation 
to the full Boltzmann equations. 

The most important effect included in our analysis is the contribution of the gauge anni- 
hilation processes, which lead to a significant reduction of the resulting lepton asymmetry for 
the low Higgs triplet mass. The dynamics of such a system is analyzed in Ref. jSj for the case 
of the conventional baryogenesis with heavy Higgs bosons in the SU(5) unification scheme. 
Our analysis is extended to the lowest possible values of the Higgs triplet mass where, as will 
be shown in the following, the annihilation effect dominates over decays and inverse decays. 
Another crucial ingredient of soft leptogenesis is the suppression of the asymmetry due to 
a small difference between boson and fermion statistics at finite temperature. We find that 
this effect for low values of the triplet mass becomes subleading compared to that due to 
soft supersymmetry-breaking terms. Let us also note that the annihilation effect becomes 
irrelevant for a triplet mass higher than about 10 10 GeV [9], for which, however, the lepton 
asymmetry in the soft leptogenesis scenario is also suppressed, as it is inversely proportional 
to the triplet mass [5|. As a result, we will conclude that the required baryon asymmetry 
can be generated only at the multi-TeV range of the Higgs triplet mass, and thus the model 
can lead to distinct collider signatures through, in particular, the production and decay of 
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a doubly charged Higgs boson 



□ 



10t lllj . This opens up another interesting possibility for 



generating the neutrino masses and mixing as well as the cosmological baryon asymmetry 
at the TeV scale, which can be tested in future colliders jl^l l^ . 



one needs to introduce a 



In the supersymmetric form of the Higgs triplet model 
vector-like pair of A = (A ++ , A + , A ) and A c = (A c , A c ~, A c0 ) with hypercharge Y = 1 
and —1, allowing for the renormalizable superpotential as follows: 

W = hLLA + XxHxHxA + X 2 H 2 H 2 A C + MAA C (1) 

where hLLA contains the neutrino mass term, huuA . The soft supersymmetry breaking 
terms relevant for us are 

-C so f t = {hA L LLA + X1A1H1H1A 

+ X 2 A 2 H 2 H 2 A C + BM AA C + h.c.} 

+ m\\A\ 2 + m\ c \A c \ 2 . (2) 

Note that we have used the same capital letters to denote the superfields as well as their 
scalar components. We will consider the universal boundary condition of soft masses; Al = 
A\ = A 2 = A and = m^c = tjiq. In the limit M 3> rriQ,A, the Higgs triplet vacuum 
expectation value (A ) = X 2 (H$) 2 /M gives the neutrino mass 

m v = 2hX 2 ^. (3) 
The mass matrix of the scalar triplets is diagonalized by 

A = ^(A + + A_) 
A'=-L(A + -A_) 

where A± are the mass eigenstates with the mass-squared values, M| = M 2 + m,Q±BM, and 
the mass-squared difference, AM 2 = IBM. In terms of the mass eigenstates, the Lagrangian 
becomes 

-C = — 7= A± [h LL + h(A L ± M) LL 

X l H l H l + X l (A l ±M)H 1 H 1 (4) 
±X* 2 H 2 H 2 ± X* 2 (A* 2 ± M) H 2 H 2 ] + h.c. 



The heavy particles A± decay to the leptonic final states, LL, LL, as well as the Higgs final 
states, HiHi, HiHi and H 2 H 2 , H 2 H 2 . Thus, the out-of-equilibrium decay A± — > LL,LL 
can lead to lepton asymmetry of the universe. 

In order to discuss how to generate a lepton asymmetry in the supersymmetric triplet 
seesaw model let us first consider the general case of a charged particle X (X) decaying 
to a final state j (j) and generating tiny CP asymmetric number densities, nx — n x and 
rij — rij. The relevant Boltzmann equations in the approximation of Maxwell-Boltzmann 
distributions are 
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where K's are the number densities in unit of the entropy density s as defined by Yx = 
rix/s ~ n x /s, Y x = {n x — n x )/s and Yj = {rij — nf)/s. Here, the CP asymmetry €j in the 
decay X — > j is defined by 



T(X - j) - T(X - j) 



(6) 



In Eq. ( |5|) . K = T x /H\ with the Hubble parameter Hi = l.QQ^fglM 2 j 'mpi at the tempera- 
ture T = M, and Bj is the branching ratio of the decay X — > j. For the relativistic degrees 
of freedom in thermal equilibrium g*, we will use the Supersymmetric Standard Model value: 
g* = 228.75. 

The evolution of the X abundance is determined by the decay and inverse decay processes, 
as well as by the annihilation effect described by the diagrams of FIG. and are accounted 
for by the functions 7d and ja, respectively. Note that the triplets are charged under the 
Standard Model gauge group and thus have nontrivial gauge annihilation effect which turns 
out to be essential in determining the final lepton asymmetry. Moreover, as a consequence 
of unitarity, the relation 2Y X + • Yj = holds, so that one can drop out the equation for 
Y x , taking the replacement: 

Y T = -lYY, (7) 



Z 3 



in the last of Eqs. In our model, the heavy particle X can be either of the six charged 
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FIG. 1: Diagrams contributing to the gauge-annihilation amplitude of triplet particles. A, A c 
represent the fermionic partners of A and A c , respectively, while A indicates a gauge boson, A a 
gaugino, h a Higgs particle, h a higgsino, / a fermion and / a sfermion. 

particles; X = A^ + , Ajjj or A^_. Each of them follows the first Boltzmann equation in Eq. (J5J) 
where 7x> and 7,4 are given by 

KAz) 
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where t„, = tan(6 l H /) with 9 W the Weinberg angle, and (3{t) = y/1 — 1~ 2 . The function j D is 
the ratio of the modified Bessel functions of the first and second kind which as usual takes into 
account the decay and inverse decay effects in the Maxwell-Boltzmann limit. The function 
7a accounts for the annihilation cross-section of a triplet component X summing all the anni- 
hilation processes; XX' — > Standard Model gauge bosons/gauginos and fermions/sfermions 
where X' is some triplet component or its fermionic partner. The separate contribution of 
each diagram in FIG. ^ is detailed in the Appendix. As far as the Standard Model part is 
concerned, our result agrees with that of Ref. with one exception: the term proportional 




FIG. 2: Decay and inverse-decay amplitudes entering in the Boltzmann equations. The solid 
lines show the decay amplitude 70 in the Maxwell-Bolzmann limit as given by Eq. (JHJ), and the 
corresponding inverse-decay amplitude 7/ d ■ The dotted and dashed curves show the result of a full 
numerical evaluation of the same amplitudes for fermionic and bosonic final states, respectively. 

to tl,, due to the mixed gauge boson (W3B) final state in diagrams (c-f) of FIG.HJ is missing 
in Ref . jsl . However, this difference concerns a subdominant contribution which is expected 
to have a negligible impact on phenomenology. The decay and inverse decay amplitudes in 
the Maxwell-Boltzmann limit are plotted in FIG. [21 along with a numerical evaluation of 
the same quantities in the case of bosonic and fermionic final states, where Bose-Einstein 
and Fermi-Dirac distributions, as well as thermal masses, are included in the calculation. 
We use this latter evaluation when we solve the full Boltzmann equations for the lepton 
asymmetry numerically. The last figure shows that the Boltzmann approximation is well 
justified as expected for the region of our relevance, z > 10. 

Given 70 and ja,, we can now analyze the thermal evolution of Y x U In FIG. EJ we plot 
the quantity (Y x — Y x )/Y x , which quantifies the departure of the triplet density from its 
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FIG. 3: Fractional departure of the triplet comoving density Yx from its equilibrium value Y x q , 
as a function of z = m/T. The higher curve shows the result of a calculation where the gauge 
annihilation effect is neglected, while the lower ones show the same quantity including annihilation 
for Im(A)=\ TeV and log 10 (M/GeV) = 8,7,6,5,4,3 from left to right. All curves are evaluated in 
the Maxwell-Boltzmann approximation. 

equilibrium value. In particular, the higher line shows the result when only the processes 
of decay and inverse decay to light particles are included in the calculation. As expected, 
since K » 1, Yx follows closely the equilibrium density Y x q with a slight deviation of 
order 10 _1 . However, annihilation is indeed important in our case and cannot be neglected. 
This is shown in the same figure by the lower curves, which represent the departure of the 
triplet density from its equilibrium value when annihilation is included. The importance of 
annihilation can be understood in the following way. The inverse decay freezes out at Zf ~ 9 

5 /2 

for K = 32 as KzJ e~ Zf = 1. On the other hand, the thermal averages of the annihilation 
and decay rate can be compared by considering the following ratio 8|: 
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where ax = KH\/M. Thus, the annihilation effect becomes negligible for M •> 10 9 GeV. 
But in our case of soft leptogenesis, higher M suppresses the lepton asymmetry as e/ oc A/M, 
so there is a tension between these two effects, and lower values of M turn out to be favored. 
In FIG. EJ one can see that, due to annihilation which freezes out at z ~ 20, Yx follows 
more closely its equilibrium density Y x q compared with the previous case, with a deviation 
which is now of order 10~ 3 . In particular, this implies that the approximation 

_yeql 

is a good one, since decoupling occurs indeed at high z. Nevertheless, in our numerical 
analysis, we solve the full Boltzmann equations where the Bose-Einstein and Fermi-Dirac 
distributions as well as thermal masses, are included properly. 

To find out the cosmo logical lepton asymmetry by the decay of X = A±, one needs to 
calculate Yj with the states j = LL and LL and thus the corresponding CP asymmetry: 

r(A± -> ll, ii) - r(A ± -> LL, 11) 

L,L y ± 

Recall that one cannot rely on the the above Boltzmann equation (jSJ) for the mechanism of 

soft leptogenesis in the supersymmetric limit of M ^> mo, A, B, as the CP asymmetries in 

the bosonic and fermionic final states takes the opposite sign, = — e^, so that the total 

asymmetry in the lepton number density vanishes, Y\ = Yf, + = 0. A non- vanishing 

lepton asymmetry arises after taking into account the supersymmetry breaking effect at 
ri 

finite temperature |5j, namely the difference between the bosonic and fermionic statistics 
given by the Bose-Einstein and Fermi-Dirac distribution, respectively. Such a thermal 
supersymmetry breaking effect can be well accounted by a slight modification of the last 
Boltzmann equation of Eq. (JHJ) resulting from the extension of the usual Maxwell-Boltzmann 
approximation to the second order, as we will show below. 

The complete form of the Boltzmann equation for the CP asymmetry in the final state j 
contains 

= J dU x dU n dU n [\A^ - \Af] 

[fx(l ± f h )(l ± f n ) - Ir ljM +/*)] + ■■• (13) 

where efll's are the phase space integration factors and Aj (Aj) is the amplitude of the decay 
X — » j (X — > j). The distribution functions fj. at thermal equilibrium are = 1/ (e^ Ei ^ — 



1) or f F . — 1/ (e^/ 71 ) + 1) for the bosonic or fermionic state ?'. U sing the effective field-theory 
approach of resummed propagators for unstable particles the effective vertices of A + 
(A + ) and the states j (j) are 

^^- ^Xffl- (14) 

where AM 2 = 2BM. For A_, one takes the interchange of + <-> — and AM 2 — > —AM 2 . 
Here, fl's are the absorptive part of two point functions; 

A; 

Calculating |A,-| 2 - \Aj\ 2 oc |S^| 2 - \S 3 X \ 2 , we get for Eq. (fT3j) . 



I^| 2 -I^i 2 ](l±^)(l±/: 



J2/ 



4 • • AM 2 

--I.n W ^n- +) (AM2)2 + n2 __ . (16) 

Here, R k include the thermal propagator effect in the cutting rule |16;j and Cj are the thermal 
phase space factor of the final states. For the bosonic and fermionic states, we have 

R B = Vl-4x B (l + f Bl + fs 2 + VbJb 2 ) 



R F = (1 - 2x^1-4x^(1 - f Fl - f F2 + 2f Fl f F2 ) 



C B = Vl-4x B (l + / Bl )(l+/ B2 ) 

C F = (1 - 2x F )Vl-4x F (l - / Fl )(l - / Fa ) (17) 

where x b ,f = m BF {T) 2 /T 2 are the thermal masses of the bosons or fermions. Let us note 
in Eq. (jlfij) that the relation 

1 



^My 3 +y 3 -v?ic)(c j R k - c k Rj) 



oc 



Iv^t i i i* k* few gl+g2 g 3+g 4 
-£lm(^iy^)(e r - e r )=0 
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holds for any final states of ji j2 and intermediate states in the loop k^^. The same is 
true for the second part of Eq. (j!3|) . In fact, this is nothing but the unitarity relation 
— > j) = J2jT(X — > j) from Eq. (|13|) . Therefore, the lepton asymmetry in the 
integrand of the Boltzmann equation (|13j) is found to be 

2|/i| 2 |Ai| 2 [Im(A L )M(|A 1 | 2 -M 2 ) 

- Im^OMd^l 2 - M 2 )]C L R Hl 
+2\h\ 2 \X 2 \ 2 {[lm(A L )M(\A 2 \ 2 - M 2 ) 

+ lm(A 2 )M(\A L \ 2 - M 2 )]C L R H2 
+ \m{A L )MMlC L R k2 + Im^OMM 2 C H2 R L } . (18) 

Note that the terms proportional to \h\ A , which do not break lepton number, disappear 
because of the previous relation of C F R~ L — C~ L R F = 0. Thus, the asymmetry in Eq. (fTHJ) 
obviously contains only the mixed terms with h\\ }2 , signaling a lepton number violation. 
With the universality condition for the soft terms (A = A^ = A\ = A 2 ), we get a simple 
equation for the lepton asymmetry as follows: 



8 1 h | 2 1 A 2 1 2 Im( A) M 3 [5 BF + 6 soft ) (19) 
where 6 BF = -[Rh 2 {C l - C L ) + C L {R &2 - R H A] 

ff^ _|_ I I 2 

and 6 so ft = Rh 2 C l — , 

putting R — C — 1 in the denominator. In the limit M 3> mo, \A\, we have 

M 2 

n±± = MT ± = —(\h\ 2 + |Ai| 2 + |A 2 | 2 ) (20) 

07T 

ignoring the small thermal effect and thus putting R^ = 1. One thus finds that, the quantity 
inside the integrand of Eq. (fH?j) is proportional to 

[5 BF + 5 SO f t ] = ei [5 BF + 5 SO ft] . (21) 

Here, one has the approximation of 5 so ft = (m^ + \ A\ 2 )/M 2 as can be seen in FIG. |U One 
can also find similar expressions for the Higgs-Higgsino final states. Recall that unitarity 
relation enforces J2j e ± — 0. The supersymmetry breaking effect 5 BF at finite temperature 
can now be encoded in the Boltzmann equation with the Maxwell-Boltzmann approximation 
by considering the expansion: l/[exp(E/T) ± 1] w exp(— E/T)[l =f exp(— E/T)]. After 
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FIG. 4: The dashed curve shows the approximation to 5bf(z) in Eq. (|22j) . while the solid line 
is the result of a numerical evaluation of the same quantity, which includes the effect of thermal 
masses and of Fermi-Dirac and Bose-Einstein distributions. The dotted curve shows the result of 
a numerical calculation of the thermal average Rh 2 Ci which multiplies the soft supersymmetry 
breaking term in Eq. (|19|) . 

the phase space integration in Eq. (|13|). one obtains the simple modification of the usual 
Boltzmann equation with the insertion of the 8bf{z) function determined by 



which gives a further suppression compared to the conventional contribution with the Bessel 
Function Ki(z). The above expression, which is monotonically decreasing in z, is valid for 
z ^> 1, and is compared in FIG.EJto a numerical calculation including the effect of thermal 
masses, which cause 5bf to vanish at small z. The latter calculation of 5bf is obtained by 
numerically evaluating the thermal average of the absorption part of the two-point function 
I1± T . Concluding the above discussions, we find that the total lepton asymmetry density 




(22) 
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Yi = Yl + Y~ L follows the approximate Boltzmann equation: 

(23) 



— = 2g A zK-f D 

az 



^(zKYx - Y?) + B t <y 9 - 2^) 



where g A = 6 counts the total number of triplet components generating the lepton asym- 
metry and S(z) = 8bf(z) + S so f t . In the above equation, the number K = T±/H 1 takes the 
minimal value of K = 32 for \h\ = | A2I 3> |Ai| as we have the relation Q|; 

2\h\\X 2 \ V0-05eV; V ; 

As one goes away from the minimum value of K with \h\ 7^ | A2 1 , the quantity q in Eq. (|21|) 
gets suppressed. Furthermore, one realizes that the resulting lepton asymmetry is maximized 
in case of B L i = B H2 g 2 ^> B H g with \h\ = | A2 1 3> |Ai|, in which case the Boltzmann 
equation for the lepton asymmetry takes the simplest form of 

dY 



2g A zKj D 

az 



i^zWx-Y^-^Y 



(25) 



; 1 

Let us now note that, taking the resonance condition B = T±, one finds the maximal value 
of €1 = , which becomes order one for A ~ M ~ TeV. 

It is now easy to find the approximate solution for Y\ from Eq. (}25|) with the insertion of 
Yx given in Eq. (jllj) . Both are found to be a fairly good approximation to the numerical 
solution of the full Boltzmann equations, as expected from our previous discussions. The 
results of our numerical calculation are shown in FIG. |3] where we plot the final lepton 
asymmetry as a function of the triplet mass M for Im(v4) from 1 to 5 TeV. When \A\ ~ M, 
one needs to recover the contributions of order |A|/M which were neglected, e.g., in Eqs. ()15|) 
and ()21|) . Taking the parameter region, mo < \A\ — lm(A), we keep those contributions in 
the numerical calculations for the curves in FIG. El A remark is in order here. For M > 10 10 
GeV, the annihilation effect becomes irrelevant and the final asymmetry is determined by 
the decay and inverse decay effects, i.e., by the value of K only, which confirms the result 
of Ref. 9|. One can see this feature in FIG.O where the lepton asymmetry as a function of 
M changes slope at about M ~ 10 10 GeV. Below this value the annihilation effect sets in, 
and the final asymmetry is strongly suppressed compared to the value one would obtain by 
extrapolating the curve with the slope for M > 10 10 GeV. 

When the trilinear coupling is larger than the triplet mass, |^4|/M > 1, besides an en- 
hancement of the CP-violating term of Eq.()21|). one could expect that the additional con- 
tribution to the coupling of the triplet particles to scalar final states enhances the total 
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FIG. 5: Final lepton asymmetry produced by triplet decay as a function of M. Different curves 
refer to Jm(A) = 1, 2, 3, 4, 5 TeV from bottom to top. 

annihilation rate, increasing substantially the value of the parameter K compared to the 
amount given by Eq. (|24jl. which is obtained in the limit \A\fM <C 1. As a consequence of 
this, the consequent additional wash-out effect could in principle suppress the ensuing lepton 
asymmetry. However, this is not the case due to the fact that, as we have already discussed, 
the annihilation process freezes out later than inverse decays, so the latter play almost no 
role in the determination of the epoch when lepton asymmetry production can start. Actu- 
ally, this epoch starts when decays eventually overcome annihilations, so a higher value of 
K can slightly anticipate it, leading so to a higher asymmetry instead than a suppression, 
although this effect is quite mild. This is what we observe in the numerical calculation 
shown in FIG. where we have assumed as before \h\ = \X 2 \ ^> \Xi\ (in order to maximize 
the amount of CP violation given by Eq. (|21jl ) and maximized the CP-violating phase (i.e., 
we have assumed Re(A) = 0). 

We also remark that sphaleron interactions are kept in thermal equilibrium even after 
the electroweak phase transition and freeze out around T sp = M/z sp ~ 90 GeV, so that 
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only the lepton asymmetry produced for z < z sp can be efficiently converted into a baryon 
asymmetry 13]. As shown in FIG. E3 due to the gauge annihilation effect, the lepton 
asymmetry production is delayed until z = Zq ~ 20. For low values of M (~ a few TeV) 
one can have z sp < z , which implies a suppression in the final lepton asymmetry. This 
explains the fast rise at low values of M of all the curves in FIG0 On the other hand, the 
dips observed in the final asymmetry for M ~ 20\A\ correspond to the case when the two 
contributions 5bf and S so f t in Eq. (}2*T]) are of the same order and cancel. This dip separates 
the two regions where Sbf or S so f t dominates in the determination of the final asymmetry. 
As shown in FIG. the CP-violating contribution from the soft supersymmetry breaking 
term S so f t in Eq. (|19|) can strongly enhance the final lepton asymmetry at low values of M. 
As a result, it is evident that the required baryon asymmetry can be reached whenever A 
and M are in the multi-TeV region. 

Before concluding our work, let us remark some experimental consequences of the model 
at future colliders. As shown above, successful baryogenesis requires a TeV-scale triplet 
mass and Yukawa couplings of the same order, h ~ A2 ~ \jm v Mjv\ ~ 10~ 6 . Thus, all the 
ow-energy lepton flavor violating processes like /1 — > or /i — > 3e are highly suppressed 



11 1 . On the other hand, future accelerators have a potential to produce such Higgs trip ets, 
in particular, the peculiar doubly charged component through the Drell-Yan processes 10|. 
Then, various features of the model can be checked by observing the branching ratios of the 
triplet decay to lepton and Higgsino pairs, in particular, A — > Ijlj , H2H2, allowing also 
to study neutrino mass patterns jll |. 

In conclusion, we have investigated baryogenesis assuming the minimal supersymmetric 
Higgs triplet model as the origin of neutrino masses and mixings. This model, with only 
one pair of triplets, can provide a mechanism for soft leptogenesis employing a CP violating 
phase and a resonant behavior in the supersymmetry breaking sector. Our analysis shows 
that the original soft leptogenesis, relying on the supersymmetry breaking effect proportional 
to the small difference between boson and fermion statistics at finite temperature cannot 
produce the right amount of baryon asymmetry due to the gauge annihilation effect. In par- 
ticular, we have calculated the full gauge-annihilation cross section including all the relevant 
supersymmetric intermediate and final states, as well as coannihilations with the fermionic 
superpartners of the triplets, finding that this effect strongly suppresses the resulting lepton 
asymmetry. On the other hand, the contribution of soft supersymmetry breaking terms, 
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particularly a sizable value for the Im(A) parameter, can enhance the lepton asymmetry 
to provide successful leptogenesis if the triplet mass is in the TeV range. In this case, the 
model predictions can be tested in future colliders by searching for a very clean signal, e.g., 
from the production and decay of doubly charged Higgs bosons. 



APPENDIX 



In this appendix we give the detailed expression for the the annihilation cross section 
shown in compact form in Eq. (jlU|) . and calculated from the diagrams (a)-(s) of FIG. ^ in 
the following, masses of light particles are neglected, while we assume a common mass M 
for the triplets and their supersymmetric partners. The reduced cross section introduced in 
Eq. (fTUj) is defined as: 

<t) = ~Y.j d^e\M\\ 

in terms of the integrated squared amplitude, averaged over the initial triplet state (hence 
the factor 1/3) and summed over the coannihilating particles, given by: 
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In the last equation we have kept within squared parentheses quantities that vanish for 
(3 — * 0, and the subscripts refer to the contributing Feynman diagrams listed in Fig. ^ In 
Eq. (|2l)j) . T a and r a /2 are the SU(2)x U(l) group generators for the triplet and doublet 
representations and f a bc are the structure constants. Assuming the minimal supersymmetric 
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Standard Model particle content, the traces are given by: 

1 



lj2tr(T a ?tr(^f = ^(14 + 11/ 



wl 1 



Y,tr{T 2 a T^)=g^ + Ul + t 



wl 1 



(27) 



ab 



and J^ab fabcfabc = 6#f . 



Since annihilation decouples for z 3> 1, the integral in Eq. 
making use of the following low-temperature expansion: 



can be approximated by 



dt 



K x {2zt) 

WM 2 



t 2 (3{t) a(t) ~ — 



49 



z 

3 /98 



where b = 47 + 32ti + — t* and b 1 = --[ — + 32t 



mi 



(28) 
(29) 



2 w 2 

Although for our results we used a numerical integration of Eq.fjHJ), we have checked that 
the above approximation leads to a good fit to the full numerical calculation for z > 10, an 
interval that safely includes the range of z relevant for the present analysis. We finally notice 
that the annihilation amplitude increases sizably in the supersymmetric theory compared to 
the Standard Model case, in significant excess of the factor 0(2) suggested by a naif expec- 
tation. In fact, the value of bo in Eq. (|2*9*|) is almost 8 times larger than the Standard Model 
value b = Q + 8t 2 w + 2t^ coming from the diagrams (c)-(f). Such an enhancement is mainly 
due to a larger number of available final states for the diagrams (i) and (1), corresponding to 
the "contact term" for scalars and to triplet-striplet annihilation to gauginos, respectively. 
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